Distributed BLAST in a grid computing context

The Basic Local Alignment Search Tool (BLAST) is one of the best known sequence comparison programs available in bioinformatics. It is used to compare query sequences to a set of target sequences, with the intention of finding similar sequences in the target set. Here, we present a distributed BLAST service which operates over a set of heterogeneous Grid resources and is made available through a Globus toolkit v.3 Grid service. This work has been carried out in the context of the BRIDGES project, a UK e-Science project aimed at providing a Grid based environment for biomedical research. Input consisting of multiple query sequences is partitioned into sub-jobs on the basis of the number of idle compute nodes available and then processed on these in batches. To achieve this, we have implemented our own Java-based scheduler which distributes sub-jobs across an array of resources utilizing a variety of local job scheduling systems

RNA-seq transcriptional profiling of peripheral blood leukocytes from cattle infected with Mycobacterium bovis

Bovine tuberculosis, caused by infection with Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including gene expression microarrays and RNA-sequencing (RNA-seq), has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix(®) GeneChip(®) Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity(®) Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top-ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA

The Impact of 18 Ancestral and Horizontally-Acquired Regulatory Proteins upon the Transcriptome and sRNA Landscape of Salmonella enterica serovar Typhimurium

Article Authors Metrics Comments Media Coverage Abstract Author Summary Introduction Results and Discussion Materials and Methods Supporting Information Acknowledgments Author Contributions References Reader Comments (0) Media Coverage (0) Figures Abstract We know a great deal about the genes used by the model pathogen Salmonella enterica serovar Typhimurium to cause disease, but less about global gene regulation. New tools for studying transcripts at the single nucleotide level now offer an unparalleled opportunity to understand the bacterial transcriptome, and expression of the small RNAs (sRNA) and coding genes responsible for the establishment of infection. Here, we define the transcriptomes of 18 mutants lacking virulence-related global regulatory systems that modulate the expression of the SPI1 and SPI2 Type 3 secretion systems of S. Typhimurium strain 4/74. Using infection-relevant growth conditions, we identified a total of 1257 coding genes that are controlled by one or more regulatory system, including a sub-class of genes that reflect a new level of cross-talk between SPI1 and SPI2. We directly compared the roles played by the major transcriptional regulators in the expression of sRNAs, and discovered that the RpoS (σ38) sigma factor modulates the expression of 23% of sRNAs, many more than other regulatory systems. The impact of the RNA chaperone Hfq upon the steady state levels of 280 sRNA transcripts is described, and we found 13 sRNAs that are co-regulated with SPI1 and SPI2 virulence genes. We report the first example of an sRNA, STnc1480, that is subject to silencing by H-NS and subsequent counter-silencing by PhoP and SlyA. The data for these 18 regulatory systems is now available to the bacterial research community in a user-friendly online resource, SalComRegulon

Effect of Stress on Viral–Bacterial Synergy in Bovine Respiratory Disease: Novel Mechanisms to Regulate Inflammation

The severity of bovine respiratory infections has been linked to a variety offactors, including environmental and nutritional changes, transportation, and socialreorganization of weaned calves. Fatal respiratory infections, however, usually occurwhen a primary viral infection compromises host defences and enhances the severityof a secondary bacterial infection. This viral–bacterial synergy can occur by a numberof different mechanisms and disease challenge models have been developed to analysehost responses during these respiratory infections. A primary bovine herpesvirus-1(BHV-1) respiratory infection followed by a secondary challenge with Mannheimia haemolyticaresults in fatal bovine respiratory disease (BRD) and host responses to these two pathogens have been studied extensively. We used this disease model todemonstrate that stress significantly altered the viral–bacterial synergy resulting infatal BRD. Functional genomic analysis revealed that BHV-1 infection enhanced toll-likereceptors (TLR) expression and increased pro-inflammatory responses whichcontribute to the severity of a Mannheimia haemolytica infection. TLRs play a criticalrole in detecting bacterial infections and inducing pro-inflammatory responses. It isdifficult to understand, however, how stress-induced corticosteroids could enhancethis form of viral–bacterial synergy. Nuclear translocation of the glucocorticoidreceptor activates cell signalling pathways which inhibit both TLR signallingand pro-inflammatory responses. The apparent conundrum between stress-inducedcorticosteroids and enhanced BRD susceptibility is discussed in terms of present data and previous investigations of stress and respiratory disease

RNA-seq Brings New Insights to the Intra-Macrophage Transcriptome of Salmonella Typhimurium.

Salmonella enterica serovar Typhimurium is arguably the world's best-understood bacterial pathogen. However, crucial details about the genetic programs used by the bacterium to survive and replicate in macrophages have remained obscure because of the challenge of studying gene expression of intracellular pathogens during infection. Here, we report the use of deep sequencing (RNA-seq) to reveal the transcriptional architecture and gene activity of Salmonella during infection of murine macrophages, providing new insights into the strategies used by the pathogen to survive in a bactericidal immune cell. We characterized 3583 transcriptional start sites that are active within macrophages, and highlight 11 of these as candidates for the delivery of heterologous antigens from Salmonella vaccine strains. A majority (88%) of the 280 S. Typhimurium sRNAs were expressed inside macrophages, and SPI13 and SPI2 were the most highly expressed pathogenicity islands. We identified 31 S. Typhimurium genes that were strongly up-regulated inside macrophages but expressed at very low levels during in vitro growth. The SalComMac online resource allows the visualisation of every transcript expressed during bacterial replication within mammalian cells. This primary transcriptome of intra-macrophage S.-Typhimurium describes the transcriptional start sites and the transcripts responsible for virulence traits, and catalogues the sRNAs that may play a role in the regulation of gene expression during infection

Lower density and shorter duration of nasopharyngeal carriage by pneumococcal serotype 1 (St217) may explain its increased invasiveness over other serotypes

ABSTRACT Streptococcus pneumoniae is a frequent colonizer of the human nasopharynx and a major cause of life-threating invasive infections such as pneumonia, meningitis and sepsis. Over 1 million people die every year due to invasive pneumococcal disease (IPD), mainly in developing countries. Serotype 1 is a common cause of IPD; however, unlike other serotypes, it is rarely found in the carrier state in the nasopharynx, which is often considered a prerequisite for disease. The aim of this study was to understand this dichotomy. We used murine models of carriage and IPD to characterize the pathogenesis of African serotype 1 (sequence type 217) pneumococcal strains obtained from the Queen Elizabeth Central Hospital in Blantyre, Malawi. We found that ST217 pneumococcal strains were highly virulent in a mouse model of invasive pneumonia, but in contrast to the generally accepted assumption, can also successfully establish nasopharyngeal carriage. Interestingly, we found that cocolonizing serotypes may proliferate in the presence of serotype 1, suggesting that acquisition of serotype 1 carriage could increase the risk of developing IPD by other serotypes. RNA sequencing analysis confirmed that key virulence genes associated with inflammation and tissue invasiveness were upregulated in serotype 1. These data reveal important new insights into serotype 1 pathogenesis, with implications for carriage potential and risk of invasive disease through interactions with other cocolonizing serotypes, an often-overlooked factor in transmission and disease progression. IMPORTANCE The pneumococcus causes serious diseases such as pneumonia, sepsis, and meningitis and is a major cause of morbidity and mortality worldwide. Serotype 1 accounts for the majority of invasive pneumococcal disease cases in sub-Saharan Africa but is rarely found during nasopharyngeal carriage. Understanding the mechanisms leading to nasopharyngeal carriage and invasive disease by this serotype can help reduce its burden on health care systems worldwide. In this study, we also uncovered the potential impact of serotype 1 on disease progression of other coinfecting serotypes, which can have important implications for vaccine efficacy. Understanding the interactions between different serotypes during nasopharyngeal carriage may lead to improved intervention methods and therapies to reduce pneumococcal invasive disease levels

Effect of Stress on Viral–Bacterial Synergy in Bovine Respiratory Disease: Novel Mechanisms to Regulate Inflammation

The severity of bovine respiratory infections has been linked to a variety of factors, including environmental and nutritional changes, transportation, and social reorganization of weaned calves. Fatal respiratory infections, however, usually occur when a primary viral infection compromises host defences and enhances the severity of a secondary bacterial infection. This viral–bacterial synergy can occur by a number of different mechanisms and disease challenge models have been developed to analyse host responses during these respiratory infections. A primary bovine herpesvirus-1 (BHV-1) respiratory infection followed by a secondary challenge with Mannheimia haemolytica results in fatal bovine respiratory disease (BRD) and host responses to these two pathogens have been studied extensively. We used this disease model to demonstrate that stress significantly altered the viral–bacterial synergy resulting in fatal BRD. Functional genomic analysis revealed that BHV-1 infection enhanced toll-like receptors (TLR) expression and increased pro-inflammatory responses which contribute to the severity of a Mannheimia haemolytica infection. TLRs play a critical role in detecting bacterial infections and inducing pro-inflammatory responses. It is difficult to understand, however, how stress-induced corticosteroids could enhance this form of viral–bacterial synergy. Nuclear translocation of the glucocorticoid receptor activates cell signalling pathways which inhibit both TLR signalling and pro-inflammatory responses. The apparent conundrum between stress-induced corticosteroids and enhanced BRD susceptibility is discussed in terms of present data and previous investigations of stress and respiratory disease

EzArray: A web-based highly automated Affymetrix expression array data management and analysis system

<p>Abstract</p> <p>Background</p> <p>Though microarray experiments are very popular in life science research, managing and analyzing microarray data are still challenging tasks for many biologists. Most microarray programs require users to have sophisticated knowledge of mathematics, statistics and computer skills for usage. With accumulating microarray data deposited in public databases, easy-to-use programs to re-analyze previously published microarray data are in high demand.</p> <p>Results</p> <p>EzArray is a web-based Affymetrix expression array data management and analysis system for researchers who need to organize microarray data efficiently and get data analyzed instantly. EzArray organizes microarray data into projects that can be analyzed online with predefined or custom procedures. EzArray performs data preprocessing and detection of differentially expressed genes with statistical methods. All analysis procedures are optimized and highly automated so that even novice users with limited pre-knowledge of microarray data analysis can complete initial analysis quickly. Since all input files, analysis parameters, and executed scripts can be downloaded, EzArray provides maximum reproducibility for each analysis. In addition, EzArray integrates with Gene Expression Omnibus (GEO) and allows instantaneous re-analysis of published array data.</p> <p>Conclusion</p> <p>EzArray is a novel Affymetrix expression array data analysis and sharing system. EzArray provides easy-to-use tools for re-analyzing published microarray data and will help both novice and experienced users perform initial analysis of their microarray data from the location of data storage. We believe EzArray will be a useful system for facilities with microarray services and laboratories with multiple members involved in microarray data analysis. EzArray is freely available from <url>http://www.ezarray.com/</url>.</p